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FREE-FLIGHT INVESTIGATION OF A 
ROCKET -PROPELLED MODEL TO DETERMINE THE AERODYNAMIC 
HEATING ON A THIN, UNSWEPT, UNTAPERED, 
MULTISPAR, ALUMINUM-ALLOY WING 
AT MACH NUMBERS UP TO 2.22 
By Emily W. Stephens 

SUMMARY 


A free-flight investigation has been made to determine some effects 
of aerodynamic heating on the structural behavior of a wing at supersonic 
speeds. The test wing was a thin, unswept, untapered, multispar, 
aluminum-alloy wing having a 20-inch chord, a 20-inch exposed semispan, 
and a circular-arc airfoil section with a thickness ratio of 5 percent. 
The wing was tested on a model propelled by a two-stage rocket-propulsion 
system to a Mach number of 2.22 and a corresponding Reynolds number per 
foot of 13.2 x 10 6 . 

Reasonably good agreement was obtained between Stanton numbers 
obtained from measured temperature -time data and values obtained by the 
theory of Van Driest for flat plates having turbulent boundary layers. 
Temperature measurements made in the skin of the wing and in the inter- 
nal structures agreed well with calculated values. 

The wing was instrumented to detect any apparent fluttering motion 
in the wing, but no evidence of flutter was observed throughout the 
flight . 


INTRODUCTION 

Results have been published (ref. l) of a previous free-flight test 
in which the test wing gave no indication of flutter during flight 
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although flutter was evident in similar wings tested in the preflight 
jet of the Langley Pilotless Aircraft Research Station at Wallops Island, 
Va. (refs. 2 and 3 ). An additional model has been flown by the Langley 
Pilotless Aircraft Research Division to obtain the effects of aerody- 
namic heating on the structural behavior of a test wing at supersonic 
speeds. The earlier flight model had a chordwise rib located at the 
midspan of the wing. In an effort to induce flutter by reducing the 
wing stiffness, the present test wing was constructed with no chordwise 
rib. The test wing, a multispar aluminum-alloy wing, was mounted as 
one of four stabilizing wings on a two-stage rocket-propelled model and 
was instrumented to obtain temperature and vibration measurements. How- 
ever, no evidence of wing flutter was recorded during flight and, there- 
fore, only aerodynamic heating data are presented in this report. The 
wings were unswept and untapered, having a 20 -inch chord, a 20 -inch 
exposed semispan, and a circular-arc airfoil section with a thickness 
ratio of 5 percent. Data were recorded up to a Mach number of 2.22 and 
a corresponding Reynolds number per foot of 13.2 x 10^. 


SYMBOLS 


b 

c 


h 

h' 

h j 

k 


l 

M 


length defined in figure 10 (c), in. 

specific heat of structural material (2024-T3 aluminum alloy), 
Btu/(lb) (°F) 

specific heat of air at constant pressure, Btu/( slug) (°F) 

local heat-transfer coefficient, Btu/(sec)(sq ft)(°F) 

effective heat-transfer coefficient across riveted joint, 
Btu/(sec)(sq ft)(°F) 

joint heat-transfer coefficient, Btu/(sec)(sq ft)(°F) 


thermal conductivity of structural material, 
sec) ( sq ft) 

length defined in figure 10 (c), in. 

Mach number 




h 

P pv 


Stanton number. 
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Np r Prandtl number 

q dynamic pressure, lb/sq ft 

R Reynolds number, — ~ 

t time, sec 

T temperature, °F 

V velocity, ft/sec 

x distance from wing leading edge (measured in free-stream 

direction), in. 

y distance from wing tip (measured normal to model center line) , 

in . 

p density, slugs/cu ft 

p viscosity, slugs/ft-sec 

t thickness, ft 

cu density of structural material (2024-T3 aluminum alloy), 

lb/cu ft 

Subscripts : 

AW adiabatic wall 

s skin 

stag stagnation 

w web 

oo free-stream conditions 

TEST VEHICLE AND TECHNIQUE 
Model 

Photographs and the general arrangement of the test vehicle are 
presented in figures 1 to 5, and the geometry and dimensions of the 
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test wing are presented in figure 4, The instrumented wing was one of 
four stabilizing wings mounted symmetrically on the test vehicle. The 
wings were the same in all respects except stiffness. The test wing 
had the same design as wings tested in the preflight jet (refs. 2 and 3 ) 
which had no chordwise ribs; however, the remaining three model wings 
were stiffened by means of three chordwise ribs per wing to minimize 
the possibility of failure of the non instrumented wings. The wings 
were unswept and untapered, having a 20 -inch chord, a 20 -inch exposed 
semispan, and a circular-arc airfoil section with a thickness ratio of 
5 percent. The wings were constructed of 2024-T3 aluminum alloy and 
had 0.064-inch-thick skins, six 0 . 025 -inch-thick internal spars, a solid 
0 . 25 -inch-thick wing-tip bulkhead, and solid leading- and trailing-edge 
sections. All rivet heads were ground flush with the wing surface and 
the entire surface of the test wing was finished to a roughness of 
approximately 35 microinches. 

Reference 1 gives a detailed description of a wing previously 
flight tested as part of this investigation. The wing of reference 1 
differed from the wing described in this report in that it had greater 
stiffness which was contributed by one chordwise rib located at the 
center of the wing. 

In addition, three wings (MW-2, MW-2-(2) , and MW-2-(3)) identical 
to the present wing were tested at M 2.0 (refs. 2 and 3) in the 
preflight jet at Wallops Island, Va. 


Test -Vehicle Instrumentation 

A detailed sketch of wing instrumentation is presented in figure 5. 

Wing temperatures were measured with six No. 30 iron-constantan 
thermocouples located 4.8 inches inboard of the wing tip. Three thermo- 
couples were located in the skin midway between the spanwise spars, one 
thermocouple was located on the center line of a web located 13.6 inches 
from the wing leading edge, one thermocouple was located on the wing-chord 
plane of the solid-wedge portion of the wing leading-edge section, and 
one thermocouple was located in a corresponding position in the solid- 
wedge portion of the wing trailing-edge section. The method used to 
install the thermocouples and record temperature measurements is given 
in detail in reference 1 . 

Since several wings identical to the flight -tested wing had either 
fluttered or, in one case, fluttered and failed during tests in the 
preflight jet (refs. 2 and 3 ); the present wing was instrumented with 
three strain gages and a vibrometer (miniaturized accelerometer) . These 
flutter detectors, located as shown in figure 5* were considered ade- 
quate to detect any fluttering motion apparent in the wing during flight. 
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The strain gages used were uncalibrated Baldwin-Lima -Hamilton EBD-1D 
gages wherein the telemeter oscillator was adjusted for maximum sensi- 
tivity in order to detect relatively small wing deflections. The telem- 
eter range was such that the ground-station galvanometer would register 
full deflection when design bending moment was applied. The vibrometer, 
though not so sensitive to small disturbances as the strain gages, was 
included as an independent means of flutter detection and was more 
suited for the measurement of violent wing motions should they occur. 

In addition, telemeter data were obtained of longitudinal accelera- 
tion and dynamic pressure. 


Flight-Test Technique 

The model was propelled by a two-stage rocket -propulsion system. 

The first stage was comprised of two 2.8-KS-9300 Cajun rocket motors 
strapped together and fired simultaneously. A JATO, 6-KS-3000, T40 
rocket motor was employed as the sustainer rocket. The model was 
intended to maintain a constant Mach number of approximately 2.0 for 
several seconds of its trajectory. This Mach number was not maintained, 
however, since the loss of a booster fin near burnout of the first-stage 
rocket disturbed the model from its intended flight path and caused the 
model to folLow a higher altitude trajectory than intended. A Mach num- 
ber of approximately 2.07 was reached at burnout of the first -stage 
rocket ( 3. 15 seconds), after which the model coasted for 1.60 seconds 
before being accelerated to a peak Mach number of approximately 2.22 at 
burnout of the second-stage sustainer rocket (10.55 seconds). The model 
followed an essentially zero-lift trajectory throughout the flight. 

A time history of velocity was obtained by a CW Doppler radar. 

Other instrumentation included an NACA modified SCR-584 radar used to 
obtain space coordinates of the model in flight and a radiosonde launched 
immediately after the test flight and tracked by a Rawin set AN/GMD-1A 
used to determine atmospheric data and wind conditions. 

Time histories of several important flight test parameters are 
presented in figure 6. An altitude of approximately 7,000 feet was 
obtained at peak Mach number. 


Precision 

The maximum errors which exist in these test data were estimated 
to be as follows: 
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. . ±5 

. . ±5 

. ±4.0 

± 0.0003 j; 

. ± 0.01 i 

i 

v 

RESULTS AND DISCUSSION t 


Temperature measurements were made with thermocouples at several 
locations in the internal structure of the wing as well as in the skin. 
These thermocouple data are presented as temperature -time histories. 
Comparisons are shown between measured temperatures and calculated tem- 
peratures. With the exception of minor wing disturbances which occurred 
at booster-fin breakup, no evidence of wing flutter was recorded during 
the flight test by either the strain gages or vibrometer. 


Skin and Internal Temperatures 

Temperature -time histories obtained from thermocouple data are pre- 
sented in figure 7« Thermocouple measurements were made at three loca- 
tions in the outer skin of the test wing, at the center line of a span- 
wise web located at the wing 0.68 chord and on the wing-chord plane of 
the solid-wedge sections of the wing leading and trailing edges. Fig- 
ure 8 shows the chordwise variations of the measured skin temperatures 
and of theoretical temperatures calculated by use of numerical integra- 
tion at several typical times. For purposes of this figure the tempera- 
ture variation near the spars has been ignored since reference 1 has 
shown the effect of conduction on surface temperatures to be small for 
the conditions of this test. The theoretical values, obtained by the 
turbulent flat-plate theory of Van Driest (ref. 4), show reasonably 
good agreement with the measured values. 

Calculations were made to estimate analytically the temperatures 
in the solid leading- and trailing -edge sections of the wing and in a 
spanwise web, and these calculated values were then compared with meas- 
ured temperatures at the same wing locations (fig. 9). The method used 
to calculate the temperatures is described in detail in the appendix. 

For purposes of calculation the cross sections of the wing were divided 
into segments as shown in figure 10 and heat-balance equations were set 
up for each block. Temperatures calculated by the digital computer for 
the internal structures of the test wing showed good agreement with 
measured values. 


T , °F . . . 
T ro , °F . . . 
V ro , ft/sec . 

p, slug/cu ft 
M 
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Heat Transfer 


The measured temperature -time data were reduced to Stanton numbers 
and are compared in figure 11 with values obtained by the theory of 
Van Driest for a flat plate with laminar and turbulent boundary layers. 
The laminar values were computed by the method of reference 5 . The 
turbulent values were computed by the method of reference 4 in which 
the Von Karman similarity law for mixing length and a Reynolds analogy 
factor based upon laminar and turbulent Prandtl numbers equal to 0.71 
and 0.86, respectively, are assumed. Theoretical values of recovery 

factor equal to ^Np r j and a constant ratio of Stanton number to 

skin-friction coefficient equal to 0.60 were used in the calculations 
of the turbulent values of Stanton number. 


Local aerodynamic conditions, obtained from measured free -stream 
data by two-dimensional shock-expansion theory, were used for both the 
experimental and theoretical calculations. No attempt was made to cor- 
rect for three-dimensional effects as previous experience has shown 
these effects to be negligible for test conditions similar to these. 


The equation used to reduce the measured temperature data to Stanton 
number is: 


AF q If 

C0 )T / (^Tau - Tc 


N 


At 


St 


AW 


pVc T 


Little credence can be attached to values of Stanton numbers 
obtained from measured temperature -time curves between 11 and 13 seconds. 
In this region the values of the forcing function T^y - T g and slope 

AT S /At are small and account for large inaccuracies in the Stanton num- 
bers. The overall trend of the measured values, though somewhat lower, 
agrees reasonably well with the theoretical turbulent -boundary-layer 
values . 


Within the accuracy of measuring heat transfer, the agreement of 
the data of the current model and the data of a similar model (ref. l) 
with turbulent theory is about the same, with slightly closer approxima- 
tion obtained for the model of reference 1. 


Comparison With Preflight-Jet Tests 

Temperature data and Stanton numbers for the flight model are com- 
pared with data obtained by the Langley Structures Research Division in 
ground tests of identical wings tested in the preflight jet of the 
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Langley Pilotless Aircraft Research Division. Jet conditions were held 
approximately constant during testing at a stagnation temperature of 
500° F and a dynamic pressure of approximately 6,000 pounds per square 
foot . 


Three wings, identical to the present test wing, were tested in 
the preflight jet. One of the ground -tested wings fluttered and failed 
during testing; a second wing fluttered near the end of the test during 
each of three runs; the third wing gave no indication of flutter. Fig- 
ure 13 presents a comparison of temperatures recorded in flight with 
temperatures recorded in the preflight jet at approximately the same 
location on the test wing. An envelope is employed to include the range 
of temperatures encountered in the preflight jet. Little temperature 
rise occurred during the first 2 seconds of the flight test when the 
model was flying at subsonic Mach numbers, whereas instantaneous heat 
rise occurred during testing in the preflight jet. Consequently, the 
temperature -time curve from the preflight jet was shifted 2 seconds to 
illustrate better the relation between the two curves during the heating 
cycle. The maximum temperatures measured in the preflight jet averaged 
100° higher than the maximum temperature attained by the flight model. 

The temperature range at which the wings fluttered in the preflight jet 
has been indicated in figure 12. These temperatures were not attained 
in flight. The heating rate for the flight model during the first half 
of sustainer firing time was less than that for the ground tests at 
comparable times. Later, the heating rates for flight and ground tests 
were comparable. 

As a matter of interest the flight and ground test data are compared 
in figure 13 using the turbulent -flow correlating factor Ng^ For 

constant Mach number, turbulent theory gives Ng^ as essentially 

constant when plotted as a function of Reynolds number. The data of this 
figure show that the heat -transfer coefficients obtained for both the 
flight and preflight- jet tests are essentially constant and of the same 
magnitude. These experimental values are lower than theoretical values 
based on the nomograph of reference 4. 


CONCLUSIONS 


Temperature -time measurements were made on a multispar, untapered, 
unswept, aluminum-alloy wing in free flight up to a Mach number of 2.22 
and a corresponding Reynolds number per foot of 13.2 X 10^. These data 
were compared with theory and with other data obtained in flight and in 
ground tests. The following conclusions were indicated: 
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1. Stanton numbers obtained from measured temperature -time data 
agree fairly well with theoretical values calculated by the theory of 
Van Driest for flat plates having turbulent boundary layers. 

2. Temperatures calculated for both the wing skin and internal 
structures closely approximate measured values. 

3. No evidence of flutter was indicated in the present test. How- 
ever, no direct comparison can be made between wings tested in the pre- 
flight jet where flutter did occur and wings tested in flight since a 
higher heating rate, higher temperatures, and a higher dynamic pressure 
were obtained in the preflight jet. 


Langley Research Center, 

National Aeronautics and 
Langley Field, Va., 


Space Administration, 
October 1, 1958- 
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APPENDIX 

CALCULATIONS OF TEMPERATURES FOR INTERNAL WING STRUCTURES 


In an effort to estimate temperatures which would result from 
flight conditions imposed in the current test, calculations have been 
made by the IBM type 704 electronic data processing machine of the 
Langley Analytical Computing Branch of the internal temperatures occurring 
at the wing-chord plane in the solid wing leading- and trailing-edge 
sections and at the center line of a spanwise web. For purposes of 
calculation the cross sections of the wing were divided as shown in 
figure 10 and heat-balance equations were set up for each block. Finite 
differences in temperature for each block were solved for by making 
simultaneous solutions of these heat-balance equations over the desired 
time range. At any given time of the calculations the material proper- 
ties were assumed to remain constant. 

The heat-balance equations included terms for heat transfer by con- 
duction and convection. The effects of radiation were not included in 
these calculations as they were considered to be negligible over the 
temperature range encountered in the flight test. 

Values of heat-transfer coefficients and adiabatic wall temperatures 
used to calculate temperatures in the wing leading- and trailing-edge 
sections were based on local flow conditions and the theory of Van Driest 
for turbulent boundary layers. In order to calculate the temperatures 
at the web center line, heat-transfer coefficients and adiabatic wall 
temperatures based on measured temperature -time data were used to deter- 
mine heat transmitted to the outer skin. Heat-transfer coefficients 
were calculated for the end blocks of the outer surface of the wing 
cross sections, and the values for intermediate blocks were obtained by 
linear interpolation. An interface conductance value (hj = 300) corre- 
sponding to a riveted aluminum-aluminum structure at the average wing 
temperature encountered during flight (T = 180° F) was chosen from 
figure 11 of reference 6. 

The following examples of heat -balance equations for several blocks 
of the wing web section (fig. 10(c)) are typical for unit spanwise 
distance: 
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(1) Block 3: 

Convection 


Conduction 



Conduction 


AT 


5 ~ 


At 

T -^b^cjoc-^ 




k 3 T s( T 2 - t 3) k 3 T s( Tk - t 3) 

^2-3 Z 5- k 


(The symbols b and Z are defined in figure 10(c) . Numbers used as 
subscripts refer to block numbers.) 


(2) Block 6: 


Convection 


Conduction 


— *— — 
Block 6 


— 

\ - 


Blodk 14 

• 4 



Conduction 


Joint conduction 


at 6 = 


At 


r 6 b 6 U)C 6 


vet 1 


AW 


k 6 r s (T 5 - T 6 ) k 6 r s (T 7 - T 6 ) 

6 ) + : + : + 

l 5-6 l 6-7 


h ' b 6 ( T lU - t 6 ) 


(see ref. 7)- 

JL + ^ + 16 

h ,i k i 4 k 6 


where h ? 
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(3) Block 16 : 


Conduction 



Conduction 


At k i 6 T w ( T l 5 _ t i 6 ) + k i 6 T wi T i 7 - t i 6 ) 

T l6 b l6P c l6 z 15-16 1 16-17 


The time increment, equal to 0.5 second, and block sizes used in 
these calculations were considered to be of appropriate size to yield 
sufficient accuracy. These values may differ, however, with different 
heating rates and material properties. 
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Figure 2 .- Rocket model on launcher. L-57-5187 . 
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Figure 5.- Instrumentation of test wing of flight model. All dimen 

sions are in inches. 
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(a) Leading -edge wedge. 

Figure 7-- Variation of measured temperatures with time for various 

chordwise stations . 
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(c) Station 2. 


Figure 7*- Continued 
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igure 7. - Concluded. 





















(b) Leading-edge wedge section 
Figure 9-- Continued. 






6 


8 


10 


12 


14 


t , sec 

(c) Trailing-edge wedge section. 
Figure 9-- Concluded. 
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Figure 15 . - Comparison of Stanton number values obtained in flight and 

in preflight jet. M = 2.0. 
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